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A large number of half-Heusler compounds have been recently proposed as three-dimensional (3D) 
topological insulators (TIs) with tunable physical properties. However, no transport measurements 
associated with the topological surface states have been observed in these half-Heusler candidates due to the 
dominating contribution from bulk electrical conductance. Here we show that, by reducing the mobility of 
bulk carriers, a two-dimensional (2D) weak antilocalization (WAL) effect, one of the hallmarks of 
topological surface states, was experimentally revealed from the tilted magnetic field dependence of 
magnetoconductance in a topologically nontrivial semimetal LuPdBi. Besides the observation of a 2D WAL 
effect, a superconducting transition was revealed at ~ UK in the same bulk LuPdBi. Quantitative 
analysis within the framework of a generalized BCS theory leads to the conclusion that the 
noncentrosymmetric superconductivity of LuPdBi is fully gapped with a possibly unconventional pairing 
character. The co-existence of superconductivity and the transport signature of topological surface states in 
the same bulk alloy suggests that LuPdBi represents a very promising candidate as a topological 
superconductor. 



Topological insulators (TIs) with time-reversal-symmetry protected helical surface states have been identified 
as a new class of quantum materials for exploiting exciting physics such as Majorana fermions and quantum 
Hall effect, as well as developing potential applications in quantum computing^"^. Very recently, a number of 
ternary half-Heusler compounds of XYZ composition (X and Y are transition or rare-earth elements and Z is a 
heavy metal) have been predicted theoretically to be three-dimensional (3D) TI candidates due to their topo- 
logically nontrivial band inversion^"^. In particular, some of half-Heusler TI candidate compounds containing 
rare-earth elements with strongly correlated /electrons already exhibit many extreme properties, such as mag- 
netism^, heavy fermion behavior and superconductivity^ which give rise to the likelihood of searching 
topological superconductors^^"^^ with a fully gapped and odd-parity Cooper pairing states in this class of materi- 
als. Despite its importance, however, making observations of the signature of the surface states through the 
electrical transport measurements in these half-Heusler TI candidate compounds is challenging as a result of the 
dominating contribution from bulk electrical conductance. Consequently, the reported large linear magnetore- 
sistance and Shubnikov-de Hass quantum oscillations in half-Heusler TIs have been ascribed to the high mobility 
of bulk carriers^^'^^. Therefore, reducing the mobility of the charge carriers in the bulk is critical to revealing the 
signature of the surface states of TIs amongst the transport properties^^"^^. As the "clean" topological surface 
transport is expected to have high mobility with Dirac-cone like dispersion, when coexisted with the low mobility 
bulk conduction, it can reveal itself by quantum-featured transport properties. In this work, LuPdBi was chosen as 
a candidate for studying TIs and other properties particularly for the following reasons. On the one hand, based 
on the band structure calculations, as presented in Figure 1, we found that the highly dispersive s-like Bi 6s band 
locates well below the Fermi level E^, whereas a less dispersive band consisting mainly of a Bi 6p band and a Lu 5d 
band is located just above the (compare ref. 26 and references therein). In such a case, the electrical transport of 
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Figure 1 | Basic electronic structure of LuPdBi. (a) The electronic band 
structure of LuPdBi with spin-orbital coupling. Tg and band is 
highlighted by blue and red line. The marked region by ellipsoid around 
corresponds to a Van Hove singularity along F-X direction; (b) Total and 
partial density of states with specific electrons, (c) A schematic illustration 
of the band inversion and hybridization. 



LuPdBi is dominated by the carriers from less dispersive pd orbitals 
and, therefore, greatly degrades the mobility of the bulk carriers. On 
the other hand, the less dispersive band exhibits a Van Hove sin- 
gularity along the F-X direction across the E^. In fact, this Van Hove 
scenario was used to explain the superconductivity in Heusler com- 
pounds^^~^°. Here, in this paper, we report magneto -transport and 
specific heat measurements performed in mm- sized LuPdBi single 
crystals. Because of the low bulk mobility, it was revealed that our 
samples displayed a two-dimensional (2D) like weak antilocalization 
(WAL) effect, one of the hallmarks of topological surface states, as 
inferred from an angle-dependent magnetoconductance study and 
quantitative analysis with 2D Hikami-Larkin-Nagaoka (HLN) the- 
ory. This provides the first experimental evidence of the signature of 
the surface states in the transport properties in half-Heusler com- 
pounds, confirming that the LuPdBi half-Heusler compound forms a 
topological nontrivial semimetal. In addition, the LuPdBi displays 
bulk superconductivity with a critical temperature of ~ 1.7 K, this 
being the highest reported so far in half-Heusler compounds. The 
nature of bulk superconductivity and the transport signature of topo- 
logical surface states reported here suggest strongly that the LuPdBi 
system represents a promising candidate for realizing a topological 
superconductor associated with intriguing Majorana edge states. 



Results and discussion 

Single crystals of LuPdBi were grown by a flux method as described in 
the Methods section. Figure 2a shows a representative powder X-ray 
diffraction pattern of the ground LuPdBi crystals. All peaks are well 
indexed to a MgAgAs-type structure with a lattice constant of 6.63 A 
as depicted in the inset, indicating that the samples crystallize in a 
single phase without impurity. The Laue back-reflection pattern 
shown in Figure 2b indicates an apparently well developed (100) 
planes, further confirming the high quality of our single crystal. In 
addition, we found that all the signals collected by EDX are from 
three elements (Lu, Pd and Bi) and an elemental mapping (Figure 2b) 
obtained by using the EDX demonstrates the uniform distribution of 
the elements Lu, Pd and Bi across the whole sample. 

Figure 3a shows the temperature dependence of longitudinal res- 
istivity of a LuPdBi single crystal obtained with standard four- 
terminal geometry from 300 to 0.3 K in a zero magnetic field. The 
values of electrical resistivity fall within the range of 0.1 to 0.4 mQ 
cm, similar to that in other members of RPdBi series^ \ and can be 
ascribed to the characteristics of a semimetal or narrow band semi- 
conductor. Although the resistivity increases monotonically upon 
cooling from 300 to 2 K, the behavior of the temperature depen- 
dences varies in the different temperature ranges. This nonmetallic 
resistivity behavior was also found in the Bi2(Se,Te)3 topological 
system with large bulk resistivity^^'^^. When further decreasing the 
temperature, the resistivity suddenly drops to zero, this being indi- 
cative of a superconducting transition at 1.7 K (see lower inset of 
Figure 3a). A close analysis of the temperature dependence of the 
resistivity reveals the following: as shown in the upper inset of 
Figure 3a, in a temperature range from 300 K down to 90 K, the 
electric conductance can be understood as a parallel circuit of an 
insulating bulk and a metallic surface. The insulating behavior of 
the bulk can be described by the three-dimensional (3D) variable- 
range hopping model (VRH): (T^^(xexp[-(T/To)-'^''Y\ When the 
temperature is below 90 K, the resistivity tends to saturate and to 
deviate dramatically from the insulating behavior described by the 
3D VRH model, which indicates strongly that a metallic transport 
channel plays a dominant role at low temperatures. 

The Hall resistivity p^y was measured in the temperature range 
from 300 to 2 K with a magnetic field up to 10 T. Shown in the inset 
of Figure 3b is representative data obtained at 2 and 300 K, respect- 
ively. The Hall resistivity curve can bend at higher fields when there 
are two kinds of charge carriers whose numbers are comparable 
participate in conduction. Howver, in our case, the Hall resistivity 
is almost linear in the whole field range, implying that the bulk 
contribution is still dominated. The positive slope and linear depend- 
ence of Pxy(H) on the magnetic field indicate that one type of charge 
carrier (hole) dominates the transport properties. The carrier density 
n and Hall mobility /Hh are then deduced based on the one- carrier 
model and the temperature dependences are presented in Figure 3b. 
The deduced values of n and fall within the range of 4.9 ~ 10.0 X 
10^^/cm^ and 300 ~ 400 cmV(V s), respectively, and decrease with 
decreasing temperatures. It should be noted that the deduced carrier 
mobility in LuPdBi is significantly lower than that in other half- 
Heusler TIs^^'^^'^^, which facilitates the detection of the signature of 
surface states in the transport properties of the samples. 

To understand the origin of the conductive channel, we measured 
the magnetoresistance (MR) of LuPdBi at different temperatures. 
Remarkably, a pronounced cusp of MR was observed in the low 
magnetic field region, as shown in the normalized MR profiles from 
2 K to 200 K in Figure 4a. This observation is reminiscent of the 
weak antilocalization (WAL) effect that has widely been reported in 
Bi-Se-Te TIs^^'^^'^^'^^. As a quantum correction to classical conduc- 
tivity, WAL can result from both the strong spin-orbit coupling in 
the bulk of sample and the topological surface states of TIs. To 
further clarify the origin of the WAL, MR measurements were per- 
formed by varying the direction in which the magnetic field was 
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Figure 2 | Structural characterizations of LuPdBi single crystal, (a) Powder XRD pattern of pulverized single -crystal LuPdBi sample (blue line), with the 
red line showing its calculation curve and the bottom trace displaying the intensity difference between the experimental and calculation data. Inset: 
crystallographic structure of LuPdBi. (b) Elemental maps of Lu (red), Pd (green), Bi (yellow) in a single crystal sample scanned by EDX. The typical mm- 
sized sample used for scanning, together with the Laue image, is shown in the left panel. 



applied to the sample. The variation of magneto- conductivity, Act = 
a(B) — G(B = 0), obtained at 10 K as a function of the perpendicular 
component of the applied magnetic field, B sin^, is shown in 
Figure 4b. It is evident that the low field data (-0.3 T to 0.3 T) 
overlaps onto a universal curve, suggesting that the magneto-con- 
ductivity depends mainly on the vertical component of the applied 
magnetic field. These data indicate strongly a 2D nature of the elec- 
trical transports that are likely to originate from the surface channel. 
If the WAL is caused mainly by spin- orbit coupling in a 3D bulk 
channel, the magneto-conductivity should be independent of the tilt 
angles of the magnetic field. 

To gain a deeper understanding of WAL phenomena, a more 
quantitative analysis is necessary. The effect of localization on the 
2D magneto -conductivity can be described by the Hikami-Larkin- 
Nagaoka (HLN) theory^^. In the limit of strong spin-orbit coupling 



and low bulk mobility regime, the weak field conductance variation 
can be written as: 



Aa(B): 



)-ln( 



4eBLl 



(1) 



Here, il/(x) is the digamma function, L^^ is the phase coherence length, 
and the parameter A = oce^/ln^hy for which in the thin film case a was 
equal to —1/2 (per surface) for WAL. By applying Eq. (1) to the 
experimental data, the fitting parameters of L^p and A are obtained. 
In Figure 4c, we show a representative fitting result of the magneto - 
conductivity curve obtained at 10 K in the weak field range — 0.3 T to 
0.3 T. It is clear that the data can be fitted well using the 2D HLN 
model. The temperature dependences of the L^p and A values are 
summarized in Figure 4d. It should be noted that the value of A in 
the order of 10° Q~\ about 10^ times larger than the theoretical 2D 
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Figure 3 | Transport properties of LuPdBi single crystal, (a) Temperature-dependence of longitudinal resistivity p^x from 0.5 to 300 K. Lower inset 
highlights the superconducting transition occurring at Tc = 1.7 K in the enlarged view ranging from 0.5 to 5 K. Upper inset shows the longitudinal 
conductivity a in function of T^^^. Red line is the fitting of the 3D variable range hopping (VRH) behavior a oc exp [ — ( T/ Tq)~^'^] to the data, in the range 
of 90-300 K; deviation from the fitting at low temperature, marked by the arrow, signifies the parallel metallic conduction, (b) Temperature dependence 
of the carrier density n (blue) and Hall mobility fi (red). Inset shows the linear behavior of the Hall resistivity Pxy^s. 5 at T = 2 and 300 K, respectively. 
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Figure 4 | Weak anit-localization (WAL) effect in LuPdBi single crystal, (a) The normalized magnetoresistance (Rxx) in perpendicular magnetic fields 
(5) at temperature ranging from 2 to 200 K. It shows the WAL feature (sharp dips at around zero magnetic field) persistent until 100 K. (b) The 
normalized magneto-conductivity (AGxx) plotted in the perpendicular magnetic field component of the magnetic field, 5sin0, at 10 K. The only 
perpendicular component dependent features indicate that WAL is induced by 2D surface electrons. Inset: schematic of the measurement setup where 6 
indicates the angle between the direction of the magnetic field and the current fiow in the sample surface, (c) The change of magneto -conductivity (Aaxx) 
at 10 K (blue circle) and its fit to the 2D WAL mode (red line), (d) Temperature dependence of phase coherence length deduced from the WAL fit 
under a field limit of 0.3 T (square). The red line shows the fit according to the Eq. (2) with p' = I and p = 2. Inset shows the temperature dependence of 
the parameter A in the Eq. (1). 



values. This disparity is most likely caused by the contribution from 
the dominated 3D bulk channels, which was also observed in another 
macroscopic system of Be2Se3 crystals^^. In addition, the nearly con- 
stant value of A with varying temperature indicates fixed constitution 
of conduction channels. The existence of 2D WAL was further con- 
firmed by examining the temperature dependence of coherence 
length L(p as shown in Figure 4d. According to the WAL of TIs^^'^^, 
both the electron -electron (e-e) scattering and electron-phonon {e- 
ph) scattering are supposed to emerge in the 3D TIs. Therefore, the 
L(p as a function of temperature may be expressed as: 



1 



1 



L2(T) L2(0) 



+ AeeTP'+AepT^ 



(2) 



Here, 1^(0) represents the zero -temperature coherence length, 
AggT^ and Ag^Prepresent the contribution from the e-e and e-ph 
interaction, respectively. We fitted Eq. (2) to the values of L(p(T) 
obtained from the fitting of Eq. (1) to the magneto- conductivity 
measured at different temperatures, as shown by the red solid line 
in Figure 4d. The best fit yields A^^= 2.64 X IQ-^nm^ K, A^p= 3.05 
X 10"^/nm^ K^. It is apparent that the temperature dependence of 
L(p(T) can be well described by Eq. (2) withp' = 1 andp= 2 and Lcp(O) 
= 94 nm, which indicates clearly that a 2D magneto-transport chan- 
nel participates in the transport of our LuPdBi crystals. 

Other strong evidence for the existence of a 2D magneto-transport 
channel in LuPdBi crystals is the observation of unusual conductance 



fluctuations under high magnetic fields (See Supplementary Material 
for additional information regarding universal conductance flucta- 
tions and magnetoresistance data). As shown in Figure S2, although 
no fixed period can be identified in the fluctuation, the gradual 
increasing of the amplitude with decreasing temperature distin- 
guishes it from a random noise signal. Recently, universal conduc- 
tance fluctuations (UCF) have been reported in Bi-Se-Te TIs which 
are likely to originate from their surface state transport^^"^\ As dis- 
cussed in the supplementary, the conduction fluctuations owned a 
2D feature that relates to the surface state, but the fluctuation ampli- 
tude is still too large comparing to the UCF characterized value of e^/ 
h. Thus, we can conclude that the experimentally observed WAL in 
LuPdBi crystals is likely to result from a combination of both the 
WAL of the 2D surface channels and WAL of 3D bulk channels. 
Moreover, we have grown and measured two other half-Heusler 
LaPdBi and YPdBi single crystals. Neither any signature of WAL 
nor conductance fluctuations in LaPdBi and YPdBi^^ were found 
in these two materials (See Supplementary Material for additional 
information regarding universal conductance fluctations and mag- 
netoresistance data). This test experiment demonstrates that, in these 
gapless half-Heusler systems, reducing the mobility of the bulk car- 
riers is an effective and simple way for probing, for example, the 2D 
surface transport properties in the presence of bulk conduction. 

The low temperature resistance data as measured under different 
magnetic fields is shown in the inset of Figure 5a. It is clear that the 
crystal becomes a superconductor below 1.7 K. To our knowledge. 
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Figure 5 | Superconducting properties of LuPdBi single crystal, (a) The upper critical field 5^2, determined as the vanishing point of the superconducting 
state (inset), was plotted with the reduced transition temperature t = T/Tc. The blue solid line shows the Ginzburg-Landau (G-L) fit to the data and the 
dotted green line represents the WHH theoretical curve. Inset: temperature dependence of resistance under different applied magnetic fields, (b) 
Temperature dependent electronic specific heat q under zero field after subtracting the phonon contribution based on the Debye fit of the normal state 
(indicated by the inset). The red solid line represents the calculated Cg/T curve given by strong-coupling BCS theory with a = A(0)/kBTc = 2.6, A(0) = 
0.32 meV, = 1.45 K. Inset: Cp/Tversus Tunder 0 and 3T magnetic fields. The 3T magnetic field is larger enough to quench superconductivity, therefore 
the 3T curve represent the normal state of the material. 



1.7 K is the highest critical temperature that has been observed in the 
family of half-Heusler alloys or in the noncentrosymmetric system. 
To gain a deeper understanding of the physics associated with super- 
conductivity, we extracted the upper critical field (T) from the 
resistance curves, as shown in Figure 5a. The value of the upper 
critical field at absolute zero temperature, 3^2(0) = 2.2 T was 
obtained by fitting the generalized Ginzburg-Landau model, Bc2(T) 
= 5,2(0)(1 )/(! + f), to the extracted data, where t = T/T^. Using 
the relation Bc2 = Oo/27i(^^, where is the flux quantum, we 
obtained superconductor coherence length, ^ = 12 nm. The mean 
free path of the charge carrier i was calculated based on 
£ = hKF/ Potie^y by assuming a spherical Fermi surface with wave 
number Kp = {^n^nY'^. Taking the experimental values n = 4.9 X 
10"' m"' and po = 4.2 X 10"" Qm (Fig. 3(b)), we obtained i = 
23 nm. Thus, £ > c^, which confirms that LuPdBi is sufficiently pure 
to allow for odd-parity SC^". We also compared the experimental 
data with the universal curve for a clean spin-singlet SC with an 
orbital limited upper- critical field 5^"2'^(0) = -0.72 x T^ldBa/dT]^ 
[Werthamer-Helfand-Hohenberg (WHH) approximation]^^. At low 
temperatures, the experimental values of Bc2{T) are clearly larger 
than those predicted by the WHH model. Actually, the ^^iCT) curve 
exhibits a nearly linear dependence over the whole temperature 
range, as for YPtBi^" and ErPdBi^'. Based on the linear dependence, 
the zero temperature upper critical field was obtained Bc2 (0) = 2.6 T 
by simple extrapolation. In both cases, the value of 5c2(0)for LuPdBi 
is well below the Pauh paramagnetic limit of 5^2(0) = 1-85 T^— 3. IT, 
suggesting an orbital pair-breaking mechanism for LuPdBi. 

To further explore the superconductivity properties of LuPdBi, we 
measured the specific heat Cp{T) as a function of temperature under 
zero and 3 T magnetic field. The experimental data are presented in 
the inset of Figure 5b as Cp{T)/T ys.T_ A jump appearing at 1.7 K 
reflects the superconducting transition. The 3 T curve corresponds 
to the normal state, because the superconductivity is completely 
quenched under 3 T magnetic field. The normal state data were fitted 
to the Debye model, Cp = + Cp^ = y^T + A^T + A^T^ and the 
fitting curve is also plotted in the same curve (the green line). The 
first term in the Debye model represents the electron specific heat, 
and the latter terms are due to phonon contribution. From the fitting, 
y„ = 11.9 mJ/molK^ A3 = 1.0 mj/mol K^ and A5 = 0.36 mJ/molK' 



are obtained. The Debye temperature ©^was estimated to be 295 K. 
Subtracting the phonon contribution from the zero-field data gives 
the electronic specific heat in the SC state, as plotted in Figure 5b. 
The sharp -jump in the electron specific heat strong indicates that the 
superconductivity in LuPdBi is fully gapped and bulk in nature. 

As discussed by Kriener et al.^^ for the topological superconductor 
CuxBi2Se3^', we analyzed the cJT data in the SC state as shown in 
Figure 5b. We attempted to fit the standard BCS model to the cJT 
data shown in Figure 5b and found that the data cannot be described 
by the standard BCS model (fitting not shown here). We then 
attempted to fit the modified BCS model proposed for strong-coup- 
ling superconductors^^ to the experimental data. The total specific 
heat Ce contains two parts: c^s (superconducting state) and Cg„ (nor- 
mal state), which can be expressed as follows: 



Ces/T-- 



4N(0) 



'hWD 



kBT 

Cen/T = yn, T>Tc 



TdA^ 
'2~dT^^ 



-yresJ<Tc 

(3) 



where E= \l 8^ + A^, AiT) is the temperature dependent gap func- 
tion. Wo is the Debye frequency and y^^^ represents a constant term 
account for contribution of the non-SC part. In the modified BCS 
model, a = A(0)/A:5Tc is taken as an adjustable parameter. By fitting 
the above model to the experiment data cJT, as clearly seen in 
Figure 5b, it can be reasonably reproduced by the modified BCS 
model (red line) with A(0) = 0.32 meV and a = 2.6. The a value 
is larger than the weak-coupling BCS value of 1 .76, indicating LuPdBi 
is a strong-coupling superconductor. We should point out that the 
pairing symmetry in LuPdBi may not be the simple isotropic s wave 
as indicated by Kriener et aV^ for the topological superconductor 
CuxBi2Se3. The fitted parameter for y^g^ is 9.5 mj/mol K", thus the 
value of ys{= yn~yres) for SC is 2.4 mJ/mol K^. It may be noted 
further that the carrier density n obtained from Hall measurements 
is 4.9 X 10'7cm' while the value of y„ is 11.9 mJ/mol for LuPdBi 
single crystal. This result seems to be inconsistent with the results 
obtained in CuxBi2Se3, where n = 1.3 X 10^7cm^ and y„ = 1.95 mj/ 
mol K^. The main reason for this disparity is likely to be due to the 
difference in the effective mass of the charge carriers. In our LuPdBi 
single crystal, the effective mass of the carriers is estimated to be 
9.5 nie {nie being the mass of free electron), which is much larger than 
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that of the topological superconductor CuxBi2Se3^^. Of course, fur- 
ther experimental and theoretical work will be necessary to reveal 
more details of the large mass enhancement in LuPdBi. 

In summary, a thorough study has been undertaken of the mag- 
neto-transport and specific heat properties of half-Heusler LuPdBi 
single crystals. A clear 2D WAL effect and universal conductance 
fluctuations, both reflecting the existence of topological surface 
states, were observed experimentally in mm-sized LuPdBi crystals. 
This is, to our knowledge, the first experimental observation of the 
transport signature of surface states in the proposed half-Heusler TIs, 
which can be ascribed to the low bulk mobility. Although 2D WAL 
was observed, due to the large contribution of the 3D bulk states, the 
bulk contribution can not excluded when connected this effect to the 
high mobility topological surface states. The band gap in the bulk 
could be opened by doping or applying pressure^ Besides the obser- 
vation of a 2D WAL effect, a superconducting transition was revealed 
at Tc ~ L7 K in the bulk LuPdBi. This is probably the highest 
transition temperature in the family of half-Heusler alloys or in the 
noncentrosymmetric system. An analysis in the framework of a gen- 
eralized BCS theory leads to the conclusion that the superconduct- 
ivity of noncentrosymmetric LuPdBi is fully gapped with a possibly 
unconventional pairing character. It will be interesting to confirm 
the fully gapped and odd-parity SC using other techniques. These 
results suggest strongly that LuPdBi is a promising candidate as a 
topological superconductor for realizing Majorana fermions for 
future application in quantum computation. 

Methods 

Sample preparation. Single crystals of LuPdBi were grown by using the self- flux 
method. The constituent elements were placed in a high-purity argon atmosphere. 
Proper ratios of high purity metals of Lu (99.95%) pieces, Pd (99.99%) grains, and Bi 
(99.999%) grains were used as starting materials. LuPdBi was grounded, mixed with 
fluxed Bi powders in an atomic ratio of 1 : 10, and placed in a Tantalum crucible, 
which was loaded into a fused quartz tube. For estimation of the best ratio of the 
fluxed material to LuPdBi, we examined the mixtures with various weight ratios, 
ranging from 20 : 1 to 1 : 1. The 10:1 mixture resulted in the best solubility and crystal 
growth. The tube was sealed under Ar gas at the pressure of 10""^ Torr and then placed 
in a furnace. We also tried to estimate the best cooling rate. The best growth 
conditions at the 10:1 ratio were as follows: The mixed sample was heated from room 
temperature to 1150°C for 24 hours, and then keeped at 1150°C for 24 hours and 
finally the sample was slowly cooled to 850°C at a rate of 2°C/hour. 

Materials characterization. The composition of the single crystal samples was 
determined by energy-dispersive X-ray (EDX) spectroscopy and the structure were 
checked by X-ray diffraction (XRD) with Cu-Ka radiation. The single-crystal 
orientation was checked by a standard Laue diffraction technique. 

Transport measurements. To measure the transport properties of a single crystal ( 1 .5 
X 1.5 X 0.2 mm^), the electrical leads in the four-probe and Hall bar configurations 
were attached onto the sample using silver paste by gold wires. The transport and 
specific heat measurements were performed in the temperature range from 0.4 to 
300 K on a Quantum Design PPMS-9 system with a ^He refrigeration insert. The 
perpendicular and longitudinal magnetoresistances were measured with the magnetic 
field perpendicular or parallel to the electrical current direction. Magnetoresistance 
versus magnetic field was also measured for several temperatures. The Hall effect was 
measured by scanning the magnetic field from — lOT to + lOT, and the Hall 
coefficient was calculated from the slope of the measured Hall effect curves. 

Band structure calculations. The electronic structure calculations in this work were 
performed using full-potential linearized augmented plane-wave method, as 
implemented in the package WIEN2K. The exchange correlation of electrons was 
treated within the local spin density approximation (LSD A) including spin orbit 
coupling (SOC). Meanwhile, al7X 17X 17 /c-point grid was used in the calculations, 
equivalent to 5000 k points in the first Brillouin zone. Moreover, the muffin-tin radii 
of the atoms are 2.5 a.u, which are generated by the system automatically. The lattice 
parameters and atomic positions were taken from our experimental data. 
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